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We calculate the masses of the nucleons and deltas in twisted mass heavy baryon chiral perturba- 
tion theory. We work to quadratic order in a power counting scheme in which we treat the lattice 
spacing, a, and the quark masses, m q , to be of the same order. We give expressions for the mass 
and the mass splitting of the nucleons and deltas both in and away from the isospin limit. We give 
an argument using the chiral Lagrangian treatment that, in the strong isospin limit, the nucleons 
remain degenerate and the delta multiplet breaks into two degenerate pairs to all orders in chiral 
perturbation theory. We show that the mass splitting between the degenerate pairs of the deltas 
first appears at quadratic order in the lattice spacing. We discuss the subtleties in the effective 
chiral theory that arise from the inclusion of isospin breaking. 
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INTRODUCTION 



Twisted mass lattice QCD (tmLQCD) Q, is an alternative regularization for lattice QCD that has recently 
received considerable attention. 1 It has the potential to match the attractive features of improved staggered fermions 
(efficient simulations absence of "exceptional configurations" 0, O(a) improvement at maximal twist operator 
mixing as in the continuum 0j 0, ) while not sharing the disadvantage of needing to take roots of the determinant 
to remove unwanted degrees of freedom. Thus tmLQCD offers a promising and interesting new way to probe the 
properties and interactions of hadrons non-perturbatively from first principles. 

Due to the limitations in computational capabilities, the quark masses, m q , used in current simulations are still 
unphysically large. Thus extrapolations in m q are necessary if physical predications are to be made from lattice 
calculations. This can be done in a systematic and model independent way through the use of chiral perturbation 
theory (%PT). Since xPT is derived in the continuum, it can be employed only after the continuum limit has been 
taken, where the lattice spacing, a, is taken to zero. However, when close to the continuum, it can be extended 
to lattice QCD at non-zero a, where discretization errors arising from the finite lattice spacing are systematically 
included in a joint expansion in a and m q . For tmLQCD with mass-degenerate quarks, the resulting "twisted mass 
| chiral perturbation theory" (tmyPT) has been formulated previously |j| , building on earlier work for the 

^ ! untwisted Wilson theory |ll IH 111- 

So far, tm^PT has only been applied to the mesonic (pionic) sector. There have been studies on pion masses and 
decay constants for m q 3> aAq CD the phase structure of tmLQCD for m q ~ a 2 Aq CD 0, ITol Il4 HM Hrl 
and quantities involving pions that do not involve final state interactions up to next-to-leading order (NLO) in the 
power counting scheme where m q ~ etAq CD 18]. However, as pointed out in Ref. |18|. many of the pionic quantities 
considered are difficult to calculate in numerical simulations because they involve quark-disconnected diagrams. This 
motivated us to extend tnrvRT to the baryon sector, which has heretofore not been done, enabling us to study 
analytically baryonic quantities that do not involve quark-disconnected diagrams. Numerical studies of the baryons 
in tmLQCD are already underway, and the first results from quenched simulations studying the nucleon and delta 
spectra have been obtained recently in Ref. . 

In the baryon sector, the extension of xPT to the lattice at finite lattice spacing to 0(a) j2jj, and to 0(a 2 ) [2]|. has 
been done for a theory with untwisted Wilson fermions. We extend that work here to include the effects of "twisting" , 
i.e. our starting underlying lattice theory is now tmLQCD. Specifically, we study the parity and flavor breaking effects 
due to twisting in the masses and mass splittings of nucleons and deltas in an SU(2) chiral effective theory. The mass 
splittings are of particular interest to us as they allow one to quantify the size of the parity-flavor breaking effects in 
tmLQCD; furthermore, they present less difficulties to numerical simulations than their counterparts in the mesonic 
sector, which involve quark-disconnected diagrams. 
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1 For a recent review see Ref. 3]. 
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We consider here tmLQCD with mass non-degenerate quarks |22j, which includes an additional parameter, the mass 
splitting, e q . This allows us to consider the theory both in and away from the strong isospin limit. With simulations 
in the near future most likely able to access the region where m q ~ aAQ CD , the power counting scheme we will adopt 
is 

1 >> e 2 ~ a A QCD ~ -^L_ ~ (1) 

AQCD Aqcd 

with e 2 denoting the small dimcnsionless expansion parameters. In the following, we will work to 0(e 4 ) in this power 
counting. 

The remainder of this article is organized as follows. In Sec. Ill Al and III Bl we briefly review the definition of 
tmLQCD with mass non-degenerate quarks, and we show how the mass splitting can be included in the Symanzik 
Lagrangian and the 0{a) meson chiral Lagrangian. Higher order corrections from the meson Lagrangian are not 
needed for the baryon observables to the order we work. In Sec. Ill CI we extend the heavy baryon xPT (HBxPT) 
to include the twisting effects to 0(e 4 )- In Sec. IIIII we present the nucleon and delta masses in tm^PT in the strong 
isospin limit, including lattice discretization errors and the flavor and parity breaking induced by the twisted mass 
term. In Sec. IIVI we extend the calculation to include isospin breaking effects and discuss the subtleties that arise. 
We conclude in Sec. [3 



II. MASS NON-DEGENERATE TWISTED MASS CHIRAL PERTURBATION THEORY 

In this section, we work out the extension of the baryon chiral Lagrangian to 0(a 2 ) given in Ref. [2l| in tmLQCD. 
We start by briefly outlining the construction of the Symanzik Lagrangian in the mass non-degenerate case, which 
follows the same procedures as those in the mass-degenerate theory with minimal modifications. 



A. The effective continuum quark level Lagrangian 

The fermionic part of the Euclidean lattice action of tmLQCD with two mass non-degenerate quarks is 

S f= ^^(^)[^7M(V*+V, I )-^V*V M + mo + i 7 5n/io-T 3 eo]^(a;), (2) 

where we have written the action given in Ref. |22T | for a general twist angle (not necessarily maximal), and in the 
so-called "twisted basis" [5|- The quark (flavor) doublets ipi and 4>i are the dimensionless bare lattice fields (with "P 
standing for lattice and not indicating left-handed), and V M and V* are the usual covariant forward and backward 
dimensionless lattice derivatives, respectively. The matrices Tj are the usual Pauli matrices acting in the flavor space, 
with T3 the diagonal matrix. The bare normal mass, too, the bare twisted mass, /iq, and the bare mass splitting eo, 
are all dimensionless parameters; an implicit identity matrix in flavor space multiplies the bare mass parameter mo. 
The notation here is that both too and eo are positive such that the upper component of the quark field is the lighter 
member of the flavor doublet with a positive bare mass. 

Note that in the mass-degenerate case, twisting can be done using any of the r^, the choice of r 3 is merely for 
convenience. Given the identity 

7T..7r. 

exp(-i-T k ) r a exp(i-T fe ) = ekabU , (3) 

one can always rotate from a basis where the twist is implemented by r a , a = 1, 2, to a basis where it is implemented 
by r 3 using the vector transformation 

$ -» ^exp(?;6>r fe ) , V -> expH6>T fe )V> , k = 1, 2, 3, 9 = ±| , (4) 

where the appropriate sign for 9 is determined by the index a. However, with T3 used here to split the quark doublet 
so that the mass term is real and flavor-diagonal, it can not be used again for twisting if the fermionic determinant 
is to remain real. 2 



2 One way to see this is to note that the mass terms mo + «75T3£to — T3eo can be written as (xq — T3J/0) exp(iQ7s) exp(i/375T3), where 
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Following the program of Symanzik 23] , and the same enumeration procedure detailed in Ref . [10| , one can obtain 
the effective continuum Lagrangian at the quark level for mass non-degenerate quarks that describes the long distance 
physics of the underlying lattice theory. Its form is constrained by the symmetries of the lattice theory. To 0(e 4 ) in 
our power counting, in which we treat aAq CD ~ m q ~ e g , we find that the Pauli term is again the only dimension 
five symmetry breaking operator just as in the mass-degenerate case |Tcj (the details of this argument are provided 
in Appendix lA")) . 

£cff = £ g + m + i7 5 Ti/z - e q T S )ip + haif) icr^F^ ip + 0(a 2 ) , (5) 

where C g is the continuum gluon Lagrangian, m is the physical quark mass, defined in the usual way by 

m = Z m (m - m c )/a, (6) 

fx is the physical twisted mass 

(i = Z M fio/a = Zp X no/a, (7) 

and e q is the physical mass splitting 

e g = Z e e /a= Z^eo/a. (8) 

The factors Zp and Z$ are matching factors for the non-singlet pseudoscalar and scalar densities respectively. Note 
that the lattice symmetries forbid additive renormalization to both /io and eo |22|. The quantity rh c is the critical 
mass, aside from an O(a) shift (see Ref. Il8j and discussion below). 

Anticipating the fact that the mesons contribute to the baryon masses only through loops, and so will be of 0(e 3 ) 
or higher, we only need to have a meson chiral Lagrangian to 0(a) for the order we work; £ e g as given in Eq. (JjjJ 
is sufficient for its construction. To build the effective chiral Lagrangian for baryons to 0[e ) on the other hand, 
terms of 0(a 2 ) in Eq. (J5J are of the appropriate size to be included. However, except for the operator which breaks 
0(4) rotation symmetry, a^ipj^D^D^D^ip, the 0(a 2 ) operators do not break the continuum symmetries in a manner 
different than the terms explicitly shown in Eq. ©, and thus their explicit form is not needed. The 0(4) breaking 
term will lead to operators in the baryon chiral Lagrangian at the order we work. However, it is invariant under 
twisting and thus contributes as those in the untwisted theory pi) . 



B. The 5(7(2) meson sector 



The low energy dynamics of the theory are described by a generalized chiral Lagrangian found by matching from 
the continuum effective Lagrangian (J5J|. As usual, the chiral Lagrangian is built from the 5(7(2) matrix- valued field 
E, which transforms under the chiral group SU(2)l x SU(2)r as 

E^LY,R\ LeSU(2) L , R€SU(2)r. (9) 

The vacuum expectation value, So = (£), breaks the chiral symmetry spontaneously down to an 5(7(2) subgroup. 
The fluctuations around E correspond to the pseudo-Goldstone bosons (pions). 

From a standard spurion analysis, the chiral Lagrangian at 0(e 2 ) is (in Euclidean space 3 ) 

C x = ^Tr(^E^St) - ^Tr( X t S + E t x ) - ^Tr(A^ + £U) , (i ) 

where / is the decay constant (normalized so that f K — 132 MeV). The quantities x an d A are spurions for the quark 
masses and discretization errors respectively. At the end of the analysis they are set to the constant values 

X — * 2B (m + ifiri - e q T 3 ) = m + ifan ~ e q r 3 , A — >2Woa=a, (11) 



xo/yo = tan/3/ tana. Thus, the twisted mass term can be transformed away leaving just the normal mass term and the mass splitting 
term. However, since this involves an (7(1) axial transformation which is anomalous, an iaFF term is introduced into the action which 
we see now is complex (because of the factor of i). Thus, since the gauge action is real, this means that the fermionic action (before the 
transformation) must be complex, and so the fermionic determinant obtained from it must also be complex. This also implies that a 
theory, where both the twist and the mass splitting are implemented by rs, is o-dependent. 
3 We will work in Euclidean space throughout this article. 
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where Bq ~ O(Aqcd) and Wq ~ C(Aq CD ) are unknown dimensionfull constants, and we have defined the quantities 
to, /t and a. 

As explained in Ref. |lSj . since the Pauli term transforms exactly as the quark mass term, they can be combined 
by using the shifted spurion 

X' = X + A, (12) 

leaving the 0(e 2 ) chiral Lagrangian unchanged from its continuum form. This corresponds at the quark level to a 
redefinition of the untwisted component of the quark mass from to to 

rri =m + aWa/B . (13) 

This shift corresponds to an 0(a) correction to the critical mass, so that it becomes 

m c = Z m fh c /a - aW /B . (14) 

Since the 0(e 2 ) Lagrangian takes the continuum form, and the mass splitting term does not contribute at this 

order, the vacuum expectation value of X at this order is that which cancels out the twist in the shifted mass matrix, 
exactly as in the mass-degenerate case: 

„ m + a + iu,T-\ . / , 

(0|E|0) iO = S = m< exp(^on) , (15) 

where 

M' = V \m + a) 2 + (j 2 . (16) 

Note that M' is the leading order result for the pion mass-squared, i.e. m 2 = M' at 0(e 2 ). If we define the physical 
quark mass by 



m q = \J ml 2 + , (17) 

then it follows from HU that 

cos Wo = m'/rriq , sinc^o = £*/ m g ■ (18) 

Details of the non-perturbative determination of the twist angle and the critical mass can be found in |l8j |. and will 
not be repeated here. 

At C(e 4 ), the mass non-degenerate chiral Lagrangian for the pions retains the same form as that in the mass 
degenerate case 0, ^| , because the mass splitting does not induce any additional symmetry breaking operators in 
£ c ff. The 0{e A ) pion Lagrangian contains the usual Gasser-Leutwyler operators of C(m 2 , mp 2 , p 4 ), where m is a 
generic mass parameter that can be to, //, or e q , as well as terms of C(am, ap 2 , a 2 ) associated with the discretization 
errors. Now as we stated earlier, since the pions will enter only through loops in typical calculations of baryon 
observables, keeping the pion masses to C(e 4 ) will lead to corrections of C(e 5 ), which is beyond the order we work. 
As our concern is not in the meson sector, the 0(e 2 ) pion Lagrangian (jlOJI is thus sufficient for our purpose in this 
work. 



C. The 517(2) baryon sector 

With the effective continuum theory and the relevant part of the effective chiral theory describing t he p ions in 



hand, we now include the lowest lying spin-i and spin-| baryons into tm%PT by using HBxPT |24l l25l l26| . which 
we will refer as the twisted mass HBxPT (tmHBxPT). In SU(2) the spin-i nucleons are described by a doublet 



N=( P n ), (19) 



and the delta resonances form a flavor quartet. As they are spin-|, they are described by a Rarita-Schwinger field, 
TV , which is totally symmetric in flavor and satisfies 7 M lj, = 0. The delta fields are normalized such that 

jilll _ ^++ jill2 _ yl22 _ _L/\0 rp222 _ ^20) 

V3 V3 
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The free Lagrangian for the nucleons and deltas to 0(e 2 ) consistent with the symmetries of the lattice theory is (in 
Euclidean space) 

C x =Niv- DN - 2 a M NM l +N -2o M NN tr{M t +)-2 a w NN tr(W+) 

+ (T^iv -DTJ + A (T M r M ) + 2 7M (T^M^TJ - 2a M {T ^) ix(M^) - 2a w (T^) tr(W+) , (21) 

where the trace is taken in flavor space, and the notation, (...), denotes contractions of the flavor (tensor) indices as 
defined in e.g. Ref. [27|. The "twisted mass" spurion field is defined by 

Mf = ~ [?m%? ± ^q)H] , m§f = ^ , (22) 

with mg° being the "twisted" mass spurion for the baryons. The "Wilson" (discretization) spurion field is defined by 

1 A 2 

W± = ^ W f ± . «>Q = , (23) 

with being the Wilson spurion for the baryons. Note that we have made simplifications using the properties of 
SU(2) matrices when writing down Eq. (|21|l . When setting the spurions to their constant values, W+ is proportional 
to the identity matrix in flavor space. Thus the operators NW+N and (T^W+T^), although allowed under the 
symmetries of tmLQCD, are not independent operators with respect to AWtr(W+) and (T M T M ) tr(W+) respectively. 
This is also true of the nucleon and delta operators involving ftA f ™ in the isospin limit (but not away from it). The 
independent operators we choose to write down are those with the simplest flavor contractions, and this will be the 
case henceforth whenever we make simplifications using the properties of SU(2). 

In Eq. (|21|l . the four- vector, v^, is the heavy baryon four- velocity, and our conventional here is that v • v = 1. The 
parameter, A, is the mass splitting between the nucleons and deltas which is independent of the quark masses (often 
referred to as the nucleon-delta mass splitting in the chiral limit), and we treat A ~ m n ~ e 2 following |24l l25l 12a. |28| . 
The chiral covariant derivative, D M , acts on the nucleon and delta fields as 

(£> M iV) i = ^AT i + (V^ j % i , 
{D^) ijk = dj% k + (V^/T^ k + (V^/T^ + (V,) k k 'T» k , . (24) 
The vector and axial-vector fields are defined by 

V, = \ + frrf) , Afi = * (^t _ gt^) , ? = E (25) 

The dimensionless low energy constants (LECs), cum, 7m, and 7tm have the same numerical values as in the 
usual untwisted two-flavor HBxPT. 

As was noted in Ref. [lg , the shifting from x to x' = X + A which corresponds to the shift of the physical mass m 
to m! at the quark level does not, in general, remove the discretization [A) term, and this is seen explicitly here with 
the presence of the discretization terms. 

The Lagrangian describing the interactions of the nucleons and deltas with the pions is 

C x = 2g A NS ■ AN - 2g AA T^S • AT^ + g AN 

The tensor eyis the rank-2 analogue of the totally antisymmetric tensor e^. The vector is the covariant spin 
operator [2J,|23. The LECs in lj2U|l are the same as those in the untwisted two-flavor HBxPT- Note that the 0(s 2 ) 
free Lagrangian i|21|l and interaction Lagrangian i|26|) are the same as those given in Ref. |2lJ when the twist is 
removed, i.e. when /i = 0. With non-vanishing twist, the mass operators carry a twisted component and the vacuum 
is "twisted" from the identity to point in the direction of the twist (the flavor ri-direction here) [Isj . 

Following Ref. 0], we expand S about its vacuum expectation value, defining the physical pion fields and the 
physical £ fields by 

E = TH ph T, Z = TZ ph V{£ ph ), T = exp(kj on /2) , H ph = exp(iV2 tv ■ r/f) , T,VeSU(2), (27) 
Now, if we make the following chiral transformation (under which the effective chiral Lagrangian is invariant) 

S^LSi? 1 , £ -> L£,V^ = V£B? , N l ^V i3 N' J , Tj? k -> V ii 'V jj 'v kk 'T% j ' k ' , 

x' — ^ L X 'R j , A -> lAr? , L , R , V S SU{2) , (28) 



TfAfe^Nu + h.c. 



i, J, 



1, 2, 3. 



(26) 
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using the particular SU{2) matrices L = W = T\ we have in the transformed effective chrial Lagrangian 



S > 2-jph 

and 



£ -» €ph , A | (tphd^lh - £} ph d^ P hj , Vf, -> - (tphd^ + tl h d^ph) , (29) 



w± - wi- = ~ ± uk^u] . w q - = ^ Tt • ( 3 °) 



Note that since L = 1$ = € SE/(2), and £ = = V% h i?, 

£ = Lt SpfcJ2 = ? = (L t^ y) . (yt^fl) = tf% h R => & - £ p/l . (31) 

We see that the £ field is now the field associated with the physical pions, and the twist is transferred from 
the twisted mass (M.*^) term to the "twisted Wilson" (VV±") term, making the mass term in the HBxPT now the 
same as that in the untwisted theory. The new mass spurion, m<g, and the "twisted Wilson" spurion, Wg", now take 
constant values 

niq — > m q - e q r 3 , Wq 1 — ► a Aq CD exp(-iw n) . (32) 

We will call this the "physical pion basis" since this is the basis where the pions are physical and we will work 
in this basis from now on, unless otherwise specified. 4 A technical point we note here is that, in the isospin limit 
where twisting can be implemented by any of the three Pauli matrices, say Tfe, the physical pion basis can be found 
following the same recipe detailed above but with n in T replaced by t\- throughout. 

Rotating to the physical pion basis where the £ field is now the physical £ p /j field in all field quantities, the form of the 
interaction Lagrangian remains unchanged as given in l|26|l . while the 0(e 2 ) free heavy baryon chiral Lagrangian 1211) 
changes to 

C x =Niv- DN -2a M ~NM +N-2a M NN tv(M + )-2a w NN tr ( W l ™ ) 

+ (T M w -DT^ + A (T fjT^) + 2 7M (T^M+T^) - 2a M (T^) tr(M+) - 2W W (T^) tr(Wf ) . (33) 

Note that W?" is also proportional to the identity matrix in flavor space when set to its constant value. Thus if we 
build the free chiral Lagrangian directly in the physical pion basis, the same simplifications due to SU (2) we used in 
writing down Eq. (|21|) apply. Note that at this point, one can not yet tell whether the nucleon (TV) and the delta (T M ) 
fields are physical. This has to be determined by the theory itself. We will return to this point when calculating the 
nucleon and delta masses below. 

At C(e 4 ), there are contributions from 0(am) and 0(a 2 ) operators. The enumeration of the operators is similar to 
that set out in Ref. pi), except now the Wilson spurion field carries a twisted component. The operators appearing 
in the 0(e 4 ) chiral Lagrangian will involve two insertions off the following: Ai±, W±", and the axial current A^. 
Note that since parity combined with flavor is conserved in tmLQCD, any one insertion of A4_ or Wi"' must be 
accompanied by another insertion of A4- or Wl 1 ". Now operators with two insertions of M + or _4 M , which contribute 
to baryon masses at tree and one-loop level respectively, have the same form as those in the untwisted theory (and 
so give the same contribution). These have been written down in [2|j and will not be repeated here. Operators with 
an insertion of either a combination of v ■ A and M,- (which have the same form as in the untwisted theory), or a 
combination of v ■ A and Wi 1 ", will also not contribute to the baryon masses at 0(e 4 ). 

At O(am), there are two independent operators contributing to the masses in the nucleon sector 



c =-JL 



rcf M+ NM+N tr(W^ ) + n™ M+ NNtv(M + ) tr(W^) 



(34) 



4 As detailed in Ref. llSl . the twist angle that one determines non-perturbatively in practice, call it lj, will differ from ojq by O(a). This 
will give rise to a relative 0(a) contribution to the pion terms. But since the pions come into baryon calculations only through loops, 
the correction will be of higher order than we work. Thus to the accuracy we work, we may use either u or uq. 
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and two independent operators contributing to the masses in the delta sector 



J1 x 



tY M+ (T^M+T^tiiW? )+t 2 



(T M T M ) tr(Wf)tr(A^+) 



(35) 



where A x = 47r/. 5 Note that there are no operators involving the commutator, [7W + ,W+ m ], because it is identically 
zero. There are also operators involving A4- <X> Wi 1 " at 0(am), but these again do not contribute to the baryon masses 
at the order we work. 

Note that the scale used in Eqs. I|34|) and [i: iofl above to make the dimensions correct is not the QCD scale, Aqcd, 
but the xPT scale, A x . This follows from the naive dimensional analysis (NDA) of Ref. |30J , and from the fact that 
these operators also function as counter terms for divergences arising from the leading loop diagrams, which have 
contributions proportional to O(am). We use the same analysis to set the dimensionfull scale in all the following 
operators which contribute to the nucleon and delta masses to the order we work. 

At 0(a 2 ), there are operators that do not break the chiral symmetry arising from the bilinear operators and four- 
quark operators (see e.g. Ref. ^| for a complete listing) in the 0{a 2 ) part of £ e ff- These give rise to the tmHBxPT 
operators 



o 2 A 3 

a ^QCD 



bNN + tiTpTJ 



(36) 



There are also chiral symmetry preserving but 0(4) rotation symmetry breaking operators which arise from the 
bilinear operator of the form a 2 ip'-f fl D ll D ll D ll ip in the 0(a 2 ) part of £ e ff. These give rise to the tmHBxPT operators 



£ v = a 2 A 3 



QCD 



b v Nv^v^v^v^N + to (TuVuVpVpVuTv) + U (T^v^v^ 



(37) 



Note that these chiral symmetry preserving operators are clearly not affected by twisting (the 0(4) symmetry breaking 
operator at the quark level from which they arise involve only derivatives with no flavor structure, and {7^, 75} = 0), 
and so they have the same form and contribute to the baryon masses in the same way as in the untwisted theory. 
The chiral symmetry breaking operators at 0(a 2 ) are those with two insertions of the Wilson spurion fields. For the 
nucleons, there are two such independent operators 



1 



nf+ NNtr(WT) tr(Wf) + n[ v - NN tr(W!™ W!™) 



W_ -TF, 



Aqcd 

and for the deltas, there are three such independent operators 
1 



A 



QCD 



*7 + (T^) tr(W^) tr(W+°) + ^~ (T^) tr(Wi w Wi w ) + T k >\W t ™f\w t ™)"'T i ' i ' k 



(38) 



(39) 



In the isospin limit where the mass splitting vanishes (e q — > 0), more simplifications occur in the C(e 4 ) chiral 
Lagrangian. The nucleon operators with coefficients and n^ M+ are the same up to a numerical factor, and 

the same holds for delta operators with coefficients i^ M+ , t 2 M+ , t^ M+ , and for operators with coefficient t 1 ~ and 



Note that in the untwisted limit, the 0(e 4 ) chiral Lagrangian reduces to that given in Ref. [2]]]. In particular, with 
the twist set to zero, operators with two insertions of W!_™ will not contribute to the nucleon or the delta mass until 
0(a 2 m) ~ 0(e 6 ). But for non- vanishing twist, they contribute at 0(a 2 ). 



III. NUCLEON AND DELTA MASSES IN THE ISOSPIN LIMIT 



There have been extensive studies of the nucleon masses, and to a lesser extent, the delta masses in HByPT and 
other variants of xPT. A partial list of references of these studies includes [S3, H^, [2^, 113, E3, US, 113, llH EM H3, 
S IS |U EJ. In this work, we are concerned with corrections to the masses of the nucleons and the deltas due 
to the effect of the twisted mass parameter. We will therefore only give expressions for the mass corrections arising 



5 We will use this definition as our convention, which differs from the more standard convention, A x = 2\Z2nf, employed by other authors. 




FIG. 1: Diagrams depicting LO and NLO mass contributions to the nucleons in tmHBxPT in the physical pion basis. The 
solid, double solid and dashed lines denote nucleons, deltas and pions respectively. The solid triangle represents an insertion 
of the twisted Wilson operator as given in Eq. l-'i'-it . The solid squares denote the couplings of the baryons to the axial current 
whose form is given in Eq. 12611 . 



from the effects of lattice discretization and twisting in tmLQCD. A calculation of the nucleon and delta masses 
in the continuum in infinite volume to 0(m 2 ) can be found in Ref. p9| . The mass corrections due to finite lattice 
spacing to 0(a 2 ) in the untwisted theory with Wilson quarks can be found in Ref. [2l|, and the leading finite volume 
modifications to the nucleon mass can be found in Ref. 

In this section, we present the results of nucleon and delta masses calculated in tmHBxPT, in the isospin limit, 
where the quark doublet is mass-degenerate, and the twist is implemented by r 3 . As we discussed in Sec. Ill Al in 
the isospin limit, the content of tmLQCD is the same regardless of which Pauli matrix is used to implement the 
twist - the action for one choice is related to another by a flavor-vector rotation. This must also hold true of the 
effective chiral theory that arises from tmLQCD. Indeed, the heavy baryon Lagrangian constructed in Sec. Ill Cl with 
Ti-twisting can be rotated into that with T3-twisting by making a vector transformation, which is given by Eq. I|28() 
but with L = R = V = exp(i|r 2 ). 



A. Nucleon Masses in the Isospin Limit 



In the continuum, the mass of the nucleons in infinite volume HBxPT with two flavor-degenerate quarks are 
organized as an expansion in powers of the quark mass, which can be written as 6 

M m = M (A, A fl ) + M$ (A, A R ) + Afg /2) (A, A R ) + Afg 3 (A, A R ) + . . . (40) 

where Ni stands for either the proton (i — p) or the neutron (i = n), and is the contribution to the i th nucleon 

of 0(ro™) calculated in the continuum and infinite volume two- flavor xPT in the isospin limit pflj . The quantity, 
Mq (A, Ar), is the renormalized nucleon mass in the chiral limit; it is independent of m q and iVj, but depends the 
renormalization scale, An, and on A, the renormalized mass splitting between the nucleons and deltas which is 
independent of the quark masses. Here, we are interested in the corrections to this formula due to the effects of 
lattice discretization and twisting arising from tmLQCD. We denote these lattice corrections to the nucleon mass at 
0(e 2n ) oj 0(m") r~j 0(a n ) (factors of Aqcd needed to make the dimensions correct are implicit here) as SM^), and 
the nucleon mass in tmHBxPT is now written as 

M$ = M Q + (M " + iJfg ) + . . . (41) 

Throughout this work, we use dimensional regularization with a modified minimal subtraction scheme where we 
consistently subtract off terms proportional to 

- - 7£ + 1 + log47r . (42) 

e 

The leading correction in tm^PT comes in at tree level, arising from the twisted Wilson nucleon operator in the 
free heavy baryon Lagrangian It reads 

6M$l {uo) = -Aa w a A 2 QCB cos(w) , (43) 

where to the accuracy we work, u> can either be loq or the twist angle non-perturbatively determined. Note that this 
correction is the same for both the proton and the neutron. At leading order, the nucleon mass is automatically 0(a) 

improved, as SM^} vanishes at maximal twist, w = 7r/2. At zero twist, to = 0, it reduces to that in the untwisted 

theory jlaEl- 



Here wc use a different convention from some of the more recent nucleon mass calculations. 
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FIG. 2: Diagrams depicting mass contributions to the nucleons at next-to-next-to-leading order (NNLO) in trnHB^PT in the 
physical pion basis. The solid, double solid and dashed lines denote nucleons, deltas and pions respectively. The solid triangle 
denotes an insertion of the twisted Wilson operator as given in <55l . The solid squares denote the coupling of the baryons to 
the axial current whose form is given in Il2(il . The clear triangle denotes a tree level insertion of the operators given in Eqs. (I34B 
and ESI - E3> - 



The next contribution to the nucleon mass comes from the leading pion loop diagrams shown in Fig. ^ However, 
at the order we work, the form of the 0(m 3 / 2 ) nucleon mass contribution is unchanged from the continuum. For 
completeness we give its full expression here 



^^-^-^'^^ (44) 
where m 2 = M' = 2B a m q , is the physical pion mass-squared given in Eq. (|16[1 . and the function T is given by 



T(m,8,A R ) = {m 2 -S 2 ) 



(2) 

The corrections to M N come from both the tree level and the one-loop diagrams, as shown in Fig. [21 The twisted 
Wilson operator in the free Lagrangian (|33|l gives rise to a tadpole diagram, which produces a contribution of O(am). 
The leading Wilson spurions also contribute to O(am) when inserted inside the pion-nucleon loops, and are partly 
cancelled by wavefunction corrections. The tree level contributions come from the operators given in Eqs. (|34H and 
(|36H - H38|) . Just as in the untwisted continuum theory, these act both as the higher dimensional operators and as 
counter terms that renormalize divergences from the lower order loop contributions. For instance, coefficients n 1 + 
and n 2 M+ are renormalized to absorb divergences from the tadpole and one-loop contributions mentioned above. As 
mentioned in Sec. Ill CI the operators with these coefficients are suppressed by A x instead of Aqcd because they are 
the counterterms for the loop divergences. These coefficients are taken to be the renormalized coefficients (finite) in 
the mass calculations, and to contain the counter terms needed in our renormalization scheme. The corrections to 
M$ read 

5M$ ( w ) = 12 a w ml log (^pj cos(u) + 16 g 2 AN {a w - a w ) [j(m„, A, A R ) + m 2 ] cos( W ) 

- 2 (< M+ + 2 nf M+ ) m q cosM - a 2 Aq CD (b + M 

, (2nf- sin 2 H -4nf+ cos 2 (w)) , (46) 



t / (5 — \ — Til" ~\~ %£ \ 

V S 2 — m 2 log — S log . , „ 

1 8 + V5 2 - m 2 + is ) & \A 2 R 



f^log(£]. «4r,-, 



' fl2 ^QCD 



where the function J is given by 



J(m, S, A R ) = (m 2 28 2 ) log ( ™ ) + 28^^ log ( \ ^ % + ^ ) . (47) 

\ a r/ \o + vo — m z + is J 

Note that the C(e 4 ) corrections are again the same for both the proton and the neutron. At maximal twist, the C(e 4 ) 
corrections are given by 

5M$ (u> = it /2) = a 2 A 3 QCB (2r%~ ~b~b v ) , (48) 

while at zero twist, these reduces to the corrections given in Ref. [2l]]. We see that the nucleon masses are also 
automatically 0(a) improved at C(e 4 )- 

To the order we work, the expressions for the nucleon mass corrections in tmHBxPT given in Eqs. (|43H and l|46|l . 
together with the untwisted continuum HBxPT expressions for the nucleon masses, provide the functional form for 
the dependence of the nucleon masses on the twist and angle, to, and the quark mass, m qi which can be used to fit 
the lattice data. 
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FIG. 3: Diagrams depicting LO and NLO mass contributions to the deltas in tmHBxPT in the physical pion basis. The solid, 
double solid and dashed lines denote nucleons, deltas and pions respectively. The solid triangle is an insertion of the twisted 
Wilson operator given in Eq. (1331 . The solid squares denote the couplings of the baryons to the axial current whose form is 
given in Eq. 12(il . 



B. Delta Masses in the Isospin Limit 

Before we present our delta mass expressions, we stress that they can only be fit to tmLQCD data for sufficiently 
large quark masses such that the delta is a stableparticle. This corresponds to > 300 MeV, which is pushing the 
bounds of validity of chiral perturbation theory |4^. However, these expressions can be used to study the convergence 
of xPT for these pion masses, where the LECs can be determined. With the value of the LECs known, the mass 
calculations can be analytically continued to pion masses where the delta becomes unstable, and be used to predict 
e.g. their lifetimes. 

The delta masses in the continuum, infinite volume HBxPT with a flavor doublet of degenerate quarks have a 
similar expansion as that for the nucleons given in Eq. (|40|) . The mass expansion of the i th delta is conventionally 
written 

M Ti = M Q (A, A H ) + A + M« (A, A„) + M% 2) (A, A„) + (A, A R ) + . . . (49) 

where 

7i = A++, T 2 = A+, T 3 = A°, r 4 = A", (50) 

and Mj^ is the contribution to the i th delta of 0(m") calculated in the continuum and infinite volume two-flavor 
%PT in the isospin limit |2flj . As in the case of the nucleons, Mo (A, Ah) is the renormalized nucleon mass in the chiral 
limit, and the parameter, A, is the renormalized mass splitting between the nucleons and deltas which is independent 
of the quark masses. 7 

Both parameters, Mq and A, are flavor singlets, and are therefore renormalized in the same way. In the nucleon 
sector, all flavor singlet mass contributions that are independent of the quark mass go into renormalizing the parameter 
Mq. In the delta sector, the choice of which parameter to renormalize, Mq or A, is arbitrary. For convenience we 
choose Mq to be the renormalized nucleon mass in the chiral limit. Thus, A is simply the difference between the 
nucleon and delta masses in the chiral limit, to a given order in the chiral expansion. Both Mq and A are parameters 
which must be fit to the lattice data. For the delta masses, we use the same renormalization scheme as for the 
nucleons. 

We again denote the correction of 0(e 2n ) to the delta masses due to the effects of lattice discretization and twisting 
arising from tmLQCD as 6Mj?\ and the delta masses in tmHBxPT are written as 

My™ = Mq + A + (m£> + 5m£>) + ... (51) 
The leading mass correction arises at tree level from the twisted Wilson delta operator given in Eq. 

5M$ (u) = -Aa w a A| CD cos(w) . (52) 
Just as for the nucleons, this does not split the delta masses and vanishes at maximal twist. 

The 0{e i ) delta mass contributions are similarly given as for the nucleons. The contributing diagrams are shown 
in Figure |3J They do not cause any splitting between the deltas, and receive no discretization corrections. For 



7 There are some subtle issues involved in calculating A in HBxPT due to the fact that it is a flavor singlet, and so can modify all 
operators/LECs in the chiral Lagrangian. However, since one can not vary A, all the LECs associated with modifying A are not 
determinable. Nevertheless, one can simply fit for the physical va lue of A determined by the lattice, and use that value in all expressions 
in which it arises. For further discussion on this point, see Refs. l29t |4(1 . 
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FIG. 4: Diagrams depicting mass contributions to the deltas at NNLO in tmHBxPT in the physical pion basis. The solid, 
double solid and dashed lines denote nucleons, deltas and pions respectively. The solid triangle denotes an insertion of the 
discretization operator in Eq. 1331 . The solid squares denote the coupling of the baryons to the axial current whose form is 
given in Eq. The clear triangle denotes a tree level insertion of the operators given in Eqs. 1351 - 1371 and 1391 . 

completeness, we list the full mass expression at this order here 

, .(3/2) _ 25 ffAA 3 _ 2 SaN t( a a s /,on 

For 772 > A, the deltas are stable particles, but for < A, the deltas become unstable (as can be seen from the 
fact that the function, !F[m v , — A, Ar), picks up an imaginary component). When the deltas are unstable, one will 
not be able to use these expressions to fit the lattice data. However, these expressions can be used fit the lattice data 
for 777^ > 300 MeV, and then analytically continued to light enough pion masses for which the deltas are unstable. 

At C(e 4 ), contributions due to the effects of twisting arise from similar diagrams as in the nucleon case, and are 
shown in Fig. 01 A splitting in the delta masses first arises at this order, which comes from the operator with coefficient 
t 2 ~ given in Ij39(l . The mass corrections read 



«Anrn o . / 777 \ \ . . 2 a7 ^QCD „i \ ( \ 

cos(w) + Ag AN (a w - crw) — t| — J(™>-k, - A, A R ) cos(w) 

X \"-R/ A x 



,,,WM + ,VKAf+\ a ^QCD / \ 2a 3 

^X 



' J — r m q cos(w) + a A QCD (t + t v ) 

+ a 2 A^ CD (itY + cos 2 (w) - 2 if" sinV) + t™~ 8 Ti sin 2 (w)) , (54) 



where 



A J" 1 f°rT 4 = A++,A- 

^ = < > 1 forT,=A+,A° ' (55) 



3 



Note the appearance of the mass splitting, S^, in SMrp . We see from above that starting at C(e 4 ), the delta multiplet 
is split into two mass-degenerate pairs, with one pair containing A ++ and A - , and the other, A + and A . At maximal 

(2) 

twist, 8Mj,.' becomes 

8m£ ] (l> = tt/2) = a 2 A 3 QCD (t + t v -2 tf~ + 8 Tt ) , (56) 

while at zero twist, it reduces to that given in Ref. [2l|. Just as in the nucleon case, the delta masses to 0(e 4 ) are 
also automatically O(a) improved. 

As is the case with the nucleons, to the order we work, the expressions for the delta mass corrections in tmHBxPT 
given in Eqs. (|52|) and (|54|) . together with the untwisted continuum HBxPT expressions for the delta masses, provide 
the functional form for the dependence of the delta masses on the twist angle, uj, and the quark mass, m q , which can 
be used to fit the lattice data. 

C. Mass splittings 

Having derived the expressions for the nucleon and delta masses in tmHBxPT to order C(e 4 ) in the isospin limit, 
we now focus on the mass splittings between the nucleons and between the deltas. The mass contributions in the 

(n) (n) 

continuum, M y N ' and My. , clearly do not give rise to mass splittings for the nucleons and deltas, since they are 
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calculated with degenerate quarks. Therefore, any mass splitting can only come from the mass corrections arising 

fromtmLQCD. 

From the results of Sec. IIII inl and Sec. IIIIBl we find that to 0(e 4 ), the protons and neutrons remain degenerate, 
while the delta multiplet splits into two degenerate pairs, with A ++ and A - in one pair, and A + and A in the other. 
The splitting between the degenerate pairs in the delta multiplet is given by 

M A+ ,o - M A++ ,- = ~ tf- a 2 A 3 QCD sin 2 H = \ a 2 A 3 QCD iL . (57) 

We reiterate here that the 0(a) uncertainty inherent in the definition of the twist angle results in a correction to 
M A +,o — M A ++,~ of 0(a 3 ) ~ C(e 6 ), which is of higher order than we work. Hence to the accuracy we work, we may 
use wo or any other non-perturbatively determined twist angle for u> above. 

Just as the case of the pion mass splitting worked out in Ref. 18], this delta splitting must vanish quadratically in 
afi = am q sin(w) on general grounds, since the masses do not violate parity. One would therefore expect, naively, the 
splitting to be 0(a 2 m 2 ) ~ 0(e 8 ). But as our results show, there is, in fact, a mass dependence in the denominator 
such that the effect is 0(e 4 ). Suppose we take a -1 = 2 GeV and Aqcd = 0.5 GeV, then we would find a mass splitting 



M A +,o - M A+ +,- = 0.04 tf- GeV. (58) 

Using naive dimensional analysis, we expect ~ C(±)? and we would have a mass splitting of the delta pairs on 
the order of 50 MeV. We point out that in the recent quenched study 0, a mass splitting of on the order of 50 to 
100 MeV is found. However, our formula can only be applied to lattice data from unquenched simulations, which have 
yet to be done, since it is derived in a fully unquenched theory. Nevertheless, it is encouraging that the quenched 
results are not dramatically different from the estimate we give above with a reasonable value of Aqce>- 

Now, the degeneracies we found for the nucleons and the delta multiplet above hold not only at C(e 4 ), but in fact 
they hold to all orders in tm^PT. This can be understood by considering the lattice Wilson-Dirac operator associated 
with the action of tmLQCD given in Eq. J2J in the isospin limit with r 3 -twisting 



D WD = j £ ^(V* + V„) - V M V A< + m + , (59) 

which has the self-adjointness property [2^ 

Ti75 D WD 7 5 T! = B\ A/D . (60) 

It follows then that the propagator for the upper and lower component of the quark doublet, i/ji(x), call them S u (x, y) 
and Sd{x,y) respectively, satisfy the relations 

75 S v (x,y)~f 5 = S' d (y,x) , 75 S d (x,y) 75 = S\(y,x) . (61) 

This means that any baryon two-point correlator which is invariant under the interchange of the quark states in the 
quark doublet combined with hermitian conjugation, leading to the degeneracies mentioned above. An argument of 
this type has been given in Ref. [l9j . 

The same can also be shown in a chiral Lagrangian treatment, as must be the case. Now one of the symmetries of 
tmLQCD with two flavor-degenerate quarks and T3-twisting is the pseudo-parity transformation, Vp, where ordinary 
parity is combined with a flavor exchange |5j, 



U (x) -> U (x P ) , x P = (-x, t) 

(62) 



U k (x) ^ Ul(x p ) , fc = l,2,3 



1pl(x) -> iTi~fo1pl(xp) 

k 4>i{x) -> -ii>i{xp)"ioTi 



where t/ M are the lattice gauge fields. At the level of HBxPT, this is manifested as the invariance of the chiral 
Lagrangian under the transformations 

p(x) <-> n(x P ) , A ++ >-(x) A-' ++ (x P ), A+>°(x) <-> A°'+(x P ), ®l =l O k (x) -> ®f =1 TiO k (ajp)ri , (63) 



where for an operator in the chiral Lagrangian, Ok is any operator matrix that contracts with the flavor indices of the 
the nucleon (N) or delta (T M ) fields in the operator. If the N or the T M fields contained in an operator have a total 



13 



of 2F flavor indices, ®^ =1 C/c is the tensor product of F operator matrices which contract with the F distinct pairs 
of these flavor indices. The degeneracies in the nucleons and the delta multiplcts discussed above would then follow 
if all the operators in the chiral Lagrangian that contribute to the baryon masses have a structure that satisfies the 
condition 

®k=iOk(x) = ®k=inO k (x P )Ti , (64) 

Consider first the case for the nucleons. Since the nucleon fields are vectors in flavor space, we can take F = 1 
without loss of generality (the nucleon fields can only couple to one operator matrix) . Since the chiral Lagrangian is 
built with just Ai±, W±", A^, and V M , the operator matrix Ok can only be constructed from combinations of these 
fields. We need not consider combinations involving just Ai+ and since the flavor structure of both is trivial, 

i.e. proportional to the identity. We also need not consider mass contributions arising from pion loops, because 
they must have the same flavor structure as the tree level local counterterms used to cancel the divergences in these 
loops. Therefore, we do not have to consider operators involving A^ and V M , which give rise to mass contributions 
only through pion- nucleon interactions. This leaves us with only combinations involving Ai- and W^™ as possible 
candidates to break the degeneracy in the nucleons. As was discussed in Sec. Ill CI because of the parity-flavor 
symmetry of tmLQCD, Ok can not contain just a single Ai- or Wl™ , but must always have an even number from the 
set {A4-, Wi™}. Now any such combination would indeed have a pure tree level part, however, it is also proportional 
to the identity in flavor space. Thus there is no operator matrix, Ok, that one can construct which violates the 
condition O^ix) = Tj.Ok(xp)Ti. 

The arguments for the case of the deltas runs similar to that for the nucleons. For the same reason given in the 
nucleon case, we need not consider operator structures that involve Ai+, Wi™, A^, and V M . We need only consider 
operator structures involving an even number from the set {Ai-, Wi"'}- For the deltas, F can be three since each 
delta field has three flavor indices. But since two of the Ok in ®\ =l Ok must come from the set {Ai-, Wi™} to satisfy 
the parity-flavor symmetry of tmLQCD, we can take F to be at most two without loss of generality. Now each of 
At- and Wi™ has a tree level part that is proportional to r 3 , thus, under V F , 0\ ® O2 where Ok can be either Ai- 
or Wl™, satisfies the symmetry condition l|64|). Therefore, one can not construct operators for the deltas that break 
the degeneracy between the pairs in the delta multiplet. 

IV. NUCLEON AND DELTA MASSES AWAY FROM THE ISOSPIN LIMIT 

In this section, we present results for mass corrections due to twisting away from the isospin limit, where the quarks 
are now mass non-degenerate. To the order we work, the corrections due to the mass splitting come in only at tree 
level. For clarity, we will only point out the change arising from the quark mass splitting; we will not repeat the 
discussion on the nucleon and delta masses that are the same both in and away from the isospin limit. 

A. The flavor-diagonal basis for the mass matrix at 0(s 4 ) 

The natural choice for splitting the quark doublet is to use the real and flavor-diagonal Pauli matrix, T3, since the 
quark states one uses on the lattice correspond to the quarks in QCD in the continuum limit. But as was discussed in 
Sec. Ill Al above, twisting can not be implemented with T3 in this case (the fermionic determinant would be complex 
otherwise), and so n is used instead. This means that the quark mass matrix in C c g given in Eq. JSJ, m+in'j^Ti — e q Ts, 
can never be made flavor-diagonal through an appropriate change of basis if both the twist and the mass splitting are 
non- vanishing, because t% and T3 can not be simultaneously diagonalized. 

Since the twist is implemented by a flavor nondiagonal Pauli matrix away from the isospin limit, flavor mixings 
are induced for non-zero twist: the quark states in tmLQCD are now linear combinations of the physical quarks of 
continuum QCD. At the level of the chiral effective theory, this manifests itself in that the hadronic states described 
by the tm%PT Lagrangian are linear combinations of the continuum QCD hadronic states we observe, viz. the pions, 
nucleons, deltas, etc. 

If the effects from twisting are perturbative as compared to the isospin breaking effects, the hadronic states described 
by tmxPT will be "perturbatively close" to their corresponding continuum QCD states, i.e. the difference between 
them is small compared to the scales in the theory (see Appendix iBl for an explicit demonstration). In this case, we 
can still extract QCD observables directly from tmxPT, as the corrections will be perturbative in the small expansion 
parameter. However, if the twisting effects are on the same order as the isospin breaking effects so that the flavor 
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mixings are large, these corrections will not be perturbative. 8 Nevertheless, one can still extract information for the 
QCD observables: One can still measure the masses of these tmLQCD hadronic states in lattice simulations, and one 
can fit these to the analytic expressions for these masses calculated in tm%PT to extract the values of the LECs. The 
LECs associated with the continuum xPT contributions have the same numerical values as in tm^PT. Therefore, if 
one determines these from tmLQCD simulations, one knows the masses of the QCD hadronic states. 

At the order we work, flavor mixings are manifested in the appearance of flavor non-conserving pion-baryon vertices 
in the Feynman rules of tmHBxPT, and in that the baryon mass matrix is not flavor-diagonal. Since we work in the 
physical pion basis where the twist is carried by the Wilson spurion (now flavor non-diagonal) instead of the mass 
spurion (now flavor-diagonal), flavor mixings can only arise from operators with one or more insertions of the Wilson 
spurion field. Because of this, the flavor non-conserving pion-baryon vertices and the non-diagonal terms in the mass 
matrix must be proportional to o, the lattice spacing, and so must vanish in the continuum limit where the effects of 
the twist are fake and can be removed by a suitable chiral change of variables 0, - 9 

For the nucleons, flavor mixings induce only unphysical flavor non-conserving pion-nucleon vertices which vanish in 
the continuum limit; the nucleon matrix is still flavor diagonal at the order we work. In fact, this is true to all orders 
in tmHBxPT. The reason is the same as that given in Sec. IIII CI We need only consider the tree level part of the 
possible operator structures that one can construct from the spurion fields in tmHBxPT. Now the only spurion field 
that has a tree level part with non-diagonal flavor structure is W_ , and as we discussed above, it must be paired 
either with another Wf™ or with which renders the flavor structure of the tree level part of the combination 

trivial. Thus we may take the basis of nucleons used in the tmHBxPT Lagrangian as the physical nucleon basis. 

For the deltas, not only are there flavor non-conserving pion-delta vertices, at the order we work, the delta mass 
matrix is already flavor nondiagonal at tree level. This happens for the deltas because the tensor nature of the T M 
field allows more freedom in the way the flavor structure of the delta operator can be constructed. Thus, in order 
to have only physical tree level mass terms for the deltas, we must change to a basis where the delta mass matrix is 
diagonal, which can now only be done order by order. 

When diagonalizing the delta mass matrix, we need, in principle, to diagonalize the mass matrix that contains 
all the mass contributions from both tree and loop level to the order that one works. But we find the difference 
between diagonalizing the delta mass matrix including both tree and loop level contributions at the order we work, 
and diagonalizing that with only the tree level mass contributions, give rise to corrections only to the loop level mass 
contributions, which are higher order than we work. Thus, we will diagonalize the delta mass matrix containing just 
the tree level mass terms in our calculation for the delta masses. 

To the order we work, if the tree level mass is given by 

v A M A v A , v A = (A++ A A+ A-) , ua = (A++ A A+ A~) T , (65) 

where v A and v A are vectors of the delta basis states used in the tmHBxPT Lagrangian, and M A is the tree level 
mass matrix, the physical delta basis is given by 

v' A = S~ 1 -v A , v' A =v A -S, (66) 

where S is the matrix of eigenvectors of M A such that 

S ■ M A ■ S^ 1 — T> , (67) 

with T> the corresponding diagonal eigenvalue matrix. This implies that 

v A M A v A = {v' A ■ S- 1 ) -{S-V- S- 1 ) ■ (S ■ v' A ) = v' A Vv' A . (68) 

The full details of the diagonalization are provided in Appendix [5] In the following sections, we will work in this 
basis for calculating the delta masses. 



A qualitative guide to the size of the flavor mixings can be found in the ratios of two-point correlation functions. Define the ratio of 
QCD delta states by 

_ (A' AJ) + (A-> A*) 
ij ~ (A* A*} + (AJ AJ> ' 

Flavor mixing should be small if the off diagonal elements of Rij are small. To determine the size of the flavor mixings quantitatively, 
one has to look at the splitting in the delta multiplet. We will discuss further in the text below. 

This shows again the convenience of the pion physical basis, where all the effects of symmetry breaking in the lattice theory are 
parametrized and contained in the Wilson spurion fields, which then vanish as the symmetries are restored in the continuum limit. 
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The Feynman rules in the new basis are obtained from the same tmHBxPT Lagrangian given above in Sec. Ill Cl but 
with each of the delta flavor states now rewritten in terms of the new delta flavor states given by the defining relations 
(166(1 . Note that changing to the new delta basis induces new unphysical flavor non-conserving vertices in the delta 
interaction terms given in l|26l) . because in terms of the new basis states, flavors are mixed. However, these flavor 
mixing components are proportional to the off-diagonal elements of S, which are proportional to the lattice spacing 
as well as the twist angle (see Appendix |BJ). Thus they vanish in the limit of vanishing twist or lattice spacing, and 
so the unphysical vertices arising from them also vanish in these limits. 

We note and reiterate here that in the isospin limit, this order by order mass matrix diagonalization is unnecessary 
as one can always rotate to a basis where the twist is flavor-diagonal from the outset, and issues of flavor nonconserving 
vertices and non-diagonal mass matrices due to flavor mixings do not arise. 10 



B. The nucleon masses 



Away from the isospin limit, the first change caused by the mass splitting occurs in the continuum mass contribution 

>s 

m u = m q — e q , m d = m q + e q , m q , e q > , (69) 



Atffl , since the quark masses 



are no longer equal. 

At the order we work, the only other change due to the mass splitting appears at 0(e 4 ) in the contribution from the 
0(am) nucleon operator with coefficient n WM+ given in Eq. I|34|) . In the isospin limit, its contribution to 8M^{oj) 
is proportional to m 9 , but away from the isospin limit, it becomes 

2 n^ M+ a A x m q cos(w) — > 2 % + a A x rrii cos(o;) , (70) 

where 

\m u iori=p 

rrn = \ . (71) 

rrid tor i = n 

The corrections to the nucleon masses from the effects of lattice discretization and twisting are otherwise the same as 
those given in Eqs. I)43|) and H46fl . 

Note that the nucleon masses are automatically 0{a) improved, just as in the isospin limit. 



C. The delta masses 



Away from the isospin limit (e q ^ 0), we calculate the delta mass and mass corrections in the basis where the delta 
mass matrix is diagonal to the order we work. This diagonalization is worked out in Appendix [Bj where we have 
obtained general expressions for the new delta basis that are valid in the range from e q = to e q ~ aAQ CD . Here, 
we present the case where e q > and e q ~ m q ~ aAg CD 3> a 2 Aq CD , which is a regime that simulations in the near 
future can probe. To the order we work, we may take the new delta basis states in this regime to be 



T{ «- |Ax) = d 



A, 



IA 



IA- 



++\ 



V3B 

f U" 

AA 1 



A ) 



rpl 



T4 «-» 



I A,) = C, 



A, 



1 + 4") | A 0)_ V3B| A++) 
4A J 1 1 AA 1 ' 



V3B 
AA 1 



(72) 



where T- — A, denote the deltas in the new basis, C,- are normalization factors, and 



A = 2e q ( 7m +^ M+ ^|^cosm) , 



B = ^-a 2 A^sin 2 H- 



(73) 



In fact, as is shown in Appendix [E] if one insists on remaining in the basis where the twist in flavor non-diagonal, one would find that 
the unphysical terms arising from flavor mixings do not vanish in the continuum limit. 
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Note that A ~ 0(e 2 ) and B ~ C(e 4 ) in our power counting, so B/A ~ C(e 2 ) and the effects of the flavor mixings is 
perturbative. 

The masses of these states are comprised of the continuum expressions given in Ref. and corrections due to the 
effects of discretization and twisting. Note the continuum expressions for the delta masses here are necessarily changed 
from that in the isospin limit because m u ^ m^. The mass corrections due to the effects of lattice discretization and 
twisting, come in at tree level; the loop contributions remain unchanged from that in the isospin limit. The tree level 
mass contributions to the order we work have been worked out in Eq. (|B7(1 . in the process of diagonalization. We list 
here the full delta mass corrections to C(e 4 ), which we denote by SMi, to the mass of the delta state denoted by T/: 



5Mi(u) = — 4 aw a A-qcd cos(co>) 
aA 2 



where 



?w " <*w) — — J{m w , - A, A R ) cos(o;) 



+ 12 % m 2 log cosH + 4 g 2 AN (a 

. ,WM + aA QCD m'i , s . ,WM+ aA QCD , s 

+ 2t x — — cos(w)+4i 2 — -r 2 — m ? cos(cj) 

+ a 2 A^ CD (i + t„) + a 2 A| CD Ut? + cos 2 (w) - 2tf - sin 2 (w) + tf- 5^ sin 2 (w)) , 



i=l,...,4, (74) 



3m„ for i = 1 

2m u + rrid for i = 2 

m u + 2md for i = 3 

Srrid for i = 4 



Si 



for z 
-I for i 



1, 4 

2, 3 



(75) 



Note that (JMi(w) as given in Eq. 
respectively, but with the changes 



f73jl. is the same as the sum of SM^ and SM^J as given in Eqs. Ij52(l and Ij54(l 



(2) 



2t 



WAL 



aA QCD 



m q cos(cl') 



„ ,WM + aA QCD m't 

2*i + — ^ 5-cos(w), 



(76) 



The full expressions for the delta masses can be obtained when the continuum contributions are included. To the 
order we work, one can obtain the complete mass expression for delta denoted by T[ to 0(e 4 ) in tmHBxPT by adding 
the mass corrections, SMi(u>), to the continuum mass of the delta denoted by Tj, whose expression can be found in 
Ref. [29|. 

We stress here that one can not take the isospin limit from any of the expressions give above in this subsection. 
They have been derived for e q 7^ and with the assumption that the twisting effects are much smaller than the isospin 
breaking effects. One must use the general formulae given in Eqs. (|B6ll and l|B7|l when considering cases where these 
conditions are not true. 

Observe that away from the isospin limit, the delta masses are also automatically 0(a) improved at maximal twist 
(u> = 7r/2), as all terms proportional to a in SM^ are proportional to cos(w) as well. Hence, to the order we work, the 
contributions due to the isospin breaking are the same as that in the continuum at maximal twist. 



V. SUMMARY 



In this paper we have studied the mass spectrum of the nucleons and the deltas in tmLQCD with mass non- 
degenerate quarks using effective field theory methods. We have extended heavy baryon chiral perturbation theory 
for SU(2) to include the effects of the twisted mass, and we have done so to 0(e 4 ) in our power counting, which 
includes operators of 0(am 2 , ap 2 , a 2 ). Using the resulting tmHBxPT, we have calculated the nucleon and the delta 
masses to C(e 4 ), and we found them to be automatically 0(a) improved as expected from the properties of tmLQCD. 

Because of the twisting, the vacuum is no longer aligned with the identity in flavor space, which has non-trivial 
effects on the physical excitations (pions) of the theory. Also, depending on whether the quarks are mass degenerate 
or not, the way twisting is implemented determines what the physical baryon states are in the theory. We have 
highlighted these subtleties when doing calculations in tmHBxPT. 

In order for the pions in the theory to be physical, we have to make a particular (non-anomalous) chiral change of 
variables to undo the twisting effects. This requires the knowledge of the twisting angle, but once that is determined, 
the physical pion basis can be determined d priori. However, whether or not the nucleons and deltas are physical must 
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still be determined from the theory. In the isospin limit, both the nucleon and the delta mass matrices are diagonal, 
and so the nucleon and delta states contained in the N and fields are physical. However, away from the isospin 
limit, only the nucleon mass matrix remains diagonal. Thus, the N field can still be regarded as physical, but the 
physical deltas are now linear combinations of the flavor states contained in the field. This can be understood from 
the fact the at the quark level, the physical QCD states, the u and d quarks, are eigenstates of T3 but not of t%. So 
only in the isospin limit, where the twist can always be implemented by the flavor-diagonal Pauli matrix, T3, are the 
states contained in the quark doublet physical quarks. Away from the isospin limit, the twist can not be implemented 
by T3 anymore, and the eigenstates of the Hamiltonian of the theory are composed of linear combinations of the u 
and d quarks. 

The physical states in tm^PT are in general a mixture of those in the untwisted X-PT. The size of the mixture is 
determined by the relative sizes of the discretization effects, which are 0(a 2 ), and the isospin splitting effects, which 
are 0(e q ). In this work, we have given general expressions for the nucleon and delta masses with respect to this 
mixing of states that are valid in the range from t q — to e q ~ oAqcd ■ 

The quantities which provided the motivation for this work and turned out to be most interesting are the mass 
splittings between the nucleons and between the deltas. We found that in the isospin limit, the nucleon masses do 
not split to any order in tm^PT, while the delta multiplet splits into two degenerate pairs. This can be understood 
from the symmetries of tmLQCD at the quark level, and as we have shown, also at the level of tm%PT. The mass 
splitting between the multiplets, M^+,o — M&++,- , first arises from a tree level contribution at C(a 2 ), and it gives an 
indication of the size of the flavor breaking in tmLQCD. This splitting in the delta multiplet will be easier to calculate 
in lattice simulations than the corresponding quantity m 2 — in the meson sector |l8j , since it involves no quark 

disconnected diagrams. 

Twisted mass HByPT can also be extended to partially quenched theories (such extension of tm%PT for pions has 
recenly been done 43]). This will be useful in the near future as numerical studies of tmLQCD move from quenched 
simulations to the more realistic, if more computationally demanding, partially quenched simulations. 



We thank Stephen Sharpe for many useful discussions, and Brian Tiburzi for helpful comments. This research was 
supported by the U.S. Department of Energy Grant Nos. DE-FG02-96ER40956 (J.M.S.W.), and DE-FG03-97ER41014 
(AW.-L.). 



APPENDIX A: ABSENCE OF ADDITIONAL DIMENSION FIVE SYMMETRY BREAKING 
OPERATORS INDUCED BY THE MASS SPLITTING 

In this appendix, we show that the mass splitting does not induce any symmetry breaking terms in the effective 
continuum Lagrangian at the quark level at quadratic order. The (mass) dimension six operators in the Symanzik 
Lagrangian we can drop for the same reason given in Ref. ^3 1 since they are either of too high order (cubic or higher 
in our expansion) or they do not break the symmetries further than those of lower dimensions. For dimension five 
operators, we will show that the only allowable terms by the symmetries of the lattice theory are those that either 
vanish by the equations of motion, or can be removed by suitable 0(a) redefinitions of the parameters in Cq, the 
effective Lagrangian in the continuum limit (the lowest order effective Lagrangian). 



In the mass-degenerate case [lOj, the only dimension five operator that appears is the Pauli term. Since in the 
limit of vanishing mass splitting (the isospin limit) the mass non-degenerate theory must be the same as the mass- 
degenerate theory, any additional operators induced by the mass splitting must be proportional to the mass splitting. 
These can only be of the form 



where the notation "P\Q" means "the set P excluding the set Q". The quantities Oo and 0\ are all the possible 
independent structures with the correct dimension, which do not lead to dimension five operators vanishing by the 
equations of motion, or are not removable by redefinitions of parameters in Co. However, none of these operators are 
allowed under the symmetries of the lattice theory. Specifically, they are forbidden by charge conjugation (C) and the 
pseudo-parity transformations that combine the ordinary parity transformation (V) with a parameter sign change 
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O = r \{l}, r - {1 , r k , 75 , l5 r k } , k = 1, 2, 3 
e^O^ : 0i = {pT , mT , nT }\{pr 3 , mr 3 } , 



dim [O ] = , 
dim [d] = 1 , 



(Al) 



V = V x (n 



(A2) 
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Structure 


C 


n 


V 






T2, 75 T 2 


Tl, 7S, 75T3 


75 Tl 


T3 


Ox 


^75 x {i, n, t 3 } 




#>75T 2 






mT 2 , m7 S T 2 


"» x {n, 7s, 7573} 








ifJ,T 2 , i/i75T 2 


ifi x {n, 75, 7 5 r 3 } 


i/iT 3 





TABLE I: The structures of the dimension five operators that are non- vanishing by the equations of motion and non-removable 
by parameter redefinitions. They are classified by the symmetries that forbid them. 



or a flavor exchange or both 

V 2 F<£q =V 2 F x{e q ^-e q ), V%, (A3) 

where 

{U (x) —> U a (xp) , xp = (-x,t) 
U k (x) ^ Ul(x p ) , k=l, 2, 3 
^(a;) -> «T2 ! 37oV'(a;p) 
$(x) -> -#(xp)7 r 2 ,3 

and C/ M are the lattice link fields. Note that we have displayed the symmetries of the lattice theory |5|,|22j in the form 
which applies to the effective continuum theory. 

In Table 1, we show explicitly which symmetry forbids each of the possible structures of Co and 0\ listed in (|A1(I . 1:L 
We group the operators in columns according to the symmetry under which they are forbidden. 

The conclusion of the above discussion is that the mass splitting does not induce any additional operators that do 
not vanish by the equations of motion, or can not be removed by redefinitions of the parameters in the theory. Thus 
beyond Cq, the effective continuum Lagrangian contains only the Pauli term to the order we work, exactly as in the 
mass-degenerate case. 



(A4) 



APPENDIX B: DIAGONALIZATION OF THE DELTA MASS MATRIX 



Here we diagonalize the tree level mass matrix for the delta states. We reiterate that the difference between first 
diagonalizing the tree level mass contributions, then calculating loop effects, versus calculating the loop contributions 
then diagonalizing, is of higher order than we work. To proceed, first we list all the independent operators to C(e 4 ) 
that have tree level mass contributions, 



0(e 2 ): iT.M + T,,). (T lt T jt ) tr(.V1 + ) . (T„T„ ) tv(U7 ) 
0(am) 

0(a 2 ) 



(T^M+T^) tr(Wf ) , (T M TJ tr(Wf ) tr(M+) 



MM/ 



(T M T M ) tr(Wf) tr(W£" ) , (T M T M ) triW^W*™) , rjf (W^) 14 ' (W tw y 3 'T^' k 



The tree level delta mass matrix at the order we work, Ma, is then given by 



v^M A v A = v A 



+ C 




v 

C14 X 4 



B 
~V3 

2B) 







c 



■ \{A + 2B)+C 



B 



A- 



VA , 



C) 



Ka ?va , 



(Bl) 



(B2) 



Most of what we show can be readily inferred from |22| . What is new here is the need for V, and the use of V% 
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where the vectors v A and i>a are vectors of the (QCD) delta basis states, 



V A = (• 



A = A++ A A+ A" 



and 



K A = 



( -A -f 3 
B 1 (A-2B) 



V3 3 





The entries in Ma are given by 



A = 2e q U + ^^f^cosH 



v A = A++ A° A+ A 



B 
V3 



o \ 





MA + 2B) -4= 



B 

A 



- 1 T 



B = tf-« 2 A^ S m 2 (,) 



(B3) 



(B4) 



C = 2 m q (7M — 2 ctm) — 4 ay/a Aq CD cos(lu) + 2 m q (t^ M+ + 2t 
(t + t v ) + a 2 A 3 QCD UtY + cos 2 (w) - 2 if - sin 2 (w) 



Mi 



oA, 



QCD 



cos(w) 



Note that to the accuracy we work, w can be either w or the non-perturbatively determined twist angle. 
Except for the operators 



(B5) 



which contribute to all other operators listed in (|B1|I above have trivial flavor structure, and so contribute to 
the identity part of Hence, to diagonalize Ma, we need only diagonalize The orthogonal matrix that 

accomplishes this is 



S 



{-2A+B+2X-) 1/2 
TT2 



/ (2A-B+2X_) 1/2 
2X 1 / 2 

V3B 

2X 1 J 2 (2A-B+2X-) 1 / 2 2X 1 _ /2 (-2A+B+2X^) 1 / 2 





2X_ 
V3B 






(2A+B+2X + ) 



1/2 







{-2A-B+2X + 



\ 



,1/2 



2 A" 



V 











+ 
V3B 



%/3-B 



(B6) 



2Xl /2 (2A+B+2X + )- 1 / 2 2X 1 + /2 {-2A-B+2X+)~ 1 / 2 J 



where X± — \/ A 2 ± AB + B 2 , and each column of S is a normalized eigenvector of Ma (and hence Ka also). The 
diagonal matrix one obtains after diagonalizing A/a is then 



V = S~ 1 ■ M A ■ S = - diag 



(-A -B- 2X_ + 3C\ 
-A-B + 2X_ + 3C 
A— B - 2X+ + 3C 
\ A- B + 2X+ + 3CJ 



(B7) 



where each entry in V is an eigenvalue of Ma . 

Now if e q 0, A 7^ 0. Hence, since in our power counting A ~ 0(e 2 ) and B ~ C(e 4 ), we may expand X± in the 
ratio of B/A ~ 0(e 2 ) < 1 as 



5 B' 2 



X ± = A X ll±- + — = A[l±-- + -—+0(e») 



IB 3B 2 



2 A 8 A 2 



(B8) 



from which it follows that 



5* = 



( 1 \ 



V o o 



V3B 
AA 



V3B 
4A 
B 

4A. 



V = diag 



( - A + C \ 

A 2B i r< 
3" IT + ° 

A 2B , r< 

"S T + ° 



(B9) 
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up to corrections of C(e 4 ) for S and 0(e 6 ) for T>. 

If e q — 0, i.e. in the isospin limit, A = and X± — B. In this case, one can not find S and T> in the isospin limit 
by taking the limit A — ► in (|B9J| . since expansion in the ratio of B/A is clearly not valid. Instead, one has to go 
back to Eq. (|B6|I and Eq. I|B7I) . which in the isospin limit reduce to 

/i 0\ 

V3 1 Q Q 

C _ 2 2 u u 

^ -i 

vo o i #y 

and the eigenvalues contained in I? given in Eq. (|B10J| are the masses of the deltas at tree level in the isospin limit 
given in Eq. (|52J) and Eq. Ij54(l . Note that as discussed in the text, in the isospin limit, we need not perform any mass 
matrix diagonalization at all, since we can simply rotate from the outset to the basis where the twist is implemented 
by the diagonal T3. 

The new delta basis states are defined by 

v' A = S^ 1 ■ v A , v' A =v A -S, (Bll) 

in which the delta mass matrix is diagonal to the order we work. By writing the old (unprimed) delta basis states in 
terms of the new (primed) basis states using the defining relations given above, i.e. 

v A = S ■ v' A , v A =v' A -S- 1 , (B12) 

the Lagrangian in the new delta basis can be obtained. Note that in the case where e q ^ 0, the new basis states 
contained in v' A are "perturbatively close" to those contained in the v A , i.e. the difference is 0(e 2 ) as can be easily 
seen from Eq. (|B9|) . This is of course not true if we are in a region where B ~ A, or e q ~ q, 2 Aq Cd . 



V = diag 



f-B + C\ 

§ +C 
-B + C 

V f +cJ 



(BIO) 



R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz, Nucl. Phys. Proc. Suppl. 83, 941 (2000). 

R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [Alpha collaboration], JHEP 0108, 058 (2001). 

R. Frezzotti, "Twisted mass lattice QCD", plenary talk at Lattice 2004, hep-lat/0409138 

A. D. Kennedy, "Algorithms for lattice QCD with dynamical fermions," plenary talk at Lattice 2004, hep-lat/0409167 

R. Frezzotti and G. C. Rossi, JHEP 0408, 007 (2004). 

C. Pena, S. Sint and A. Vladikas, JHEP 0409, 069 (2004). 

R. Frezzotti and G. C. Rossi, JHEP 0410, 070 (2004). 

G. Minister and C. Schmidt, Europhys. Lett. 66, 652 (2004). 

L. Scorzato, Eur. Phys. J. C 37, 445 (2004). 

S. R. Sharpe and J. M. S. Wu, Phys. Rev. D 70, 094029 (2004). 

S. Sharpe and R. Singleton, Jr, Phys. Rev. D 58, 074501 (1998). 

G. Rupak and N. Shoresh, Phys. Rev. D 66, 054503 (2002). 

O. Bar, G. Rupak and N. Shoresh, Phys. Rev. D 70, 034508 (2004). 

G. Minister, JHEP 0409, 035 (2004). 

G. Minister, C. Schmidt and E. E. Scholz, |hep4at /0409066 
S. Sharpe and J. Wu, Nucl. Phys. Proc. Suppl. 140, 323 (2005). 
S. Aoki and O. Bar, Phys. Rev D 70, 116011 (2 004). 
S. R. Sharpe and J. M. S. Wu, |hep-lat/0411021| 

A. M. Abdel-Rehim, R. Lewis and R. M. Woloshyn, hep-lat/0503007 
S. R. Beane a nd M. J. Savage, Phys. Rev. D 68, 114502 (2003) 

B. C. Tiburzi, hep-lat/0501020 

R. Frezzotti and G. U. Rossi, Nucl. Phys. Proc. Suppl. 128, 193 (2004). 
K. Symanzik, Nucl. Phys. B 226, 187 (1983a); B 227, 205 (1983b). 

E. Jenkins and A. V. Manohar, UCSD-PTH-91-30, talk presented at the Workshop on Effective Field Theories of the 

Standard Model, Dobogoko, Hungary, Aug 1991. 

E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 

E. Jenkins, Nucl. Phys. B 368, 190 (1992). 

J. N. Labrenz and S. R. Sharpe, Phys. Rev. D 54, 4595 (1996). 

T. R. Hemmert, B. R. Holstein and J. Kambor, J. P hys. G 24, 1831 (1998). 

B. C. Tiburzi and A. Walker-Loud, hep-lat/0501018 

A. Manohar and H. Georgi, Nucl. Phys. B 234, 189 (1984). 



21 



[31] V. Bernard, N. Kaiser, and U. G. Meissner, Z. Phys. C 60, 111 (1993). 

[32] R. F. Lebed, Nucl. Phys. B 430, 295 (1994). 

[33] R. F. Lebed, Phys. Lett. B 329, 479 (1994). 

[34] M. K. Banerjee and J. Milana, Phys. Rev. D 52, 6451 (1995). 

[35] B. Borasoy and U. G. Meissner, Annals Phys. 254, 192 (1997). 

[36] T. Becher and H. Leutwyler, Eur. Phys. J. C 9, 643 (1999). 

[37] M. Frink and U. G. Meissner, JHEP 0407, 028 (2004). 

[38] A. Walker-Loud, Nucl. Phys. A 747, 476 (2005). 

[39] B. C. Tiburzi and A. Walker-Loud, Nucl. Phys. A 748, 513 (2005). 

[40] B. C. Tiburzi, Phys. Rev. D 71, 034501 (2 005). 

[41] B. C. Lehnhart, J. Gegelia, and S. Scherer, hep-ph/0412092 
[42] S. R. Beane, Nucl. Phys. B 695, 192 (2004). 

[43] G. Munster, C. Schmidt and E. E. Scholz, Europhys. Lett. 86, 639 (2004) [arXiv:hep-lat /0402003 . 



